Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  595 1715 (N48E44 L=345)   550    109        1 ( 571/1685- 612/1746) -35.5 E47-E40 N45-N50
      2  831 1497 (N32E16 L= 13)   580    172        1 ( 759/1446- 893/1543) -24.4 E21-E10 N29-N35
      3  427 1390 (N23E45 L=344)   682    138        1 ( 403/1331- 449/1467) -77.5 E47-E42 N19-N28
      4 1370 1322 (N22W22 L= 52)  1681    384        3 (1207/1231-1535/1420) -26.6 W11-W33 N16-N28
      5  893 1316 (N21E10 L= 20)   755    223        4 ( 852/1226- 941/1399) -65.6 E12-E05 N15-N26
      6 1042 1302 (N20W00 L= 29)   224     77        2 (1019/1271-1072/1339) -63.4 E01-W02 N18-N22
      7  672 1106 (N07E23 L=  7)  1014    243        1 ( 618/1003- 746/1223) -55.3 E26-E17 N00-N13
      8  811 1021 (N02E13 L= 16)   328    128        0 ( 749/ 981- 863/1069)  39.8 E16-E10 N00-N05
      9  555  946 (S04E30 L=359)   255     75        1 ( 522/ 927- 592/ 963)  21.8 E32-E27 S04-S02
     10 1302  255 (S52W32 L= 61)   422     88        0 (1280/ 233-1329/ 280) -32.0 W28-W36 S54-S48
  

    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