Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  563 1577 (N15E43 L=141)   275    113        0 ( 522/1557- 608/1601)  -7.1 E47-E39 N14-N16
      2  321 1483 (N06E60 L=124)  1500    300        9 ( 243/1411- 376/1573)  57.5 E65-E55 N00-N11
      3 1686 1293 (N26W36 L=220)   542     77        2 (1655/1272-1723/1321)   0.0 W32-W38 N23-N27
      4 1368 1119 (N06W17 L=201)   268     81        0 (1333/1087-1403/1144)  59.0 W15-W18 N03-N08
      5  859 1032 (S10E09 L=175)   231     84        0 ( 842/ 986- 875/1074) -53.1 E11-E08 S13-S08
      6  620  909 (S22E22 L=162)   634    251        1 ( 510/ 854- 731/ 973)  37.6 E29-E15 S27-S16
      7 1559  467 (S23W52 L=236)   764    220        2 (1470/ 409-1612/ 524)  68.2 W48-W54 S28-S18
      8 1411  492 (S27W40 L=224)   480    109        2 (1358/ 465-1477/ 522)  50.2 W37-W43 S29-S23
      9  865  205 (S66W24 L=208)   458     54        1 ( 836/ 185- 891/ 227)  18.4 W20-W27 S67-S62
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y