Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  846 1466 (N31E14 L= 38)   404     74        1 ( 816/1439- 891/1485)  21.8 E16-E10 N29-N32
      2 1536 1319 (N22W36 L= 87)   291     67        0 (1501/1296-1564/1341) -26.6 W33-W37 N20-N23
      3  409 1228 (N14E43 L=  8)   717    200        0 ( 380/1138- 454/1308) -51.3 E47-E38 N08-N19
      4 1129 1250 (N17W06 L= 58)   343     70        0 (1085/1231-1167/1267)  11.3 W03-W08 N16-N18
      5  407  844 (S10E42 L= 10)   313     52        2 ( 381/ 821- 430/ 881)  76.0 E44-E39 S11-S07
      6 1763  849 (S08W53 L=105)   346     39        0 (1741/ 824-1788/ 871)  21.8 W51-W56 S09-S06
      7  241  825 (S12E59 L=353)   898    101        2 ( 200/ 788- 274/ 852) -20.6 E63-E54 S14-S10
      8  738  745 (S15E18 L= 34)   223     55        0 ( 713/ 723- 759/ 773)  45.0 E19-E16 S16-S12
      9  789  577 (S26E15 L= 36)   302     73        0 ( 757/ 545- 822/ 606) -36.9 E17-E13 S28-S23
     10 1431  423 (S37W34 L= 86)  1084    211        2 (1401/ 336-1455/ 516) -65.0 W30-W39 S44-S29
  

    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