Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1260 1853 (N65W44 L=348)  1131    151        2 (1192/1829-1316/1879)   2.0 W31-W63 N62-N67
      2 1561 1577 (N37W50 L=354)   444     55        1 (1529/1562-1595/1593)   7.1 W45-W53 N36-N38
      3  967 1383 (N27E03 L=301)  1087    383        1 ( 908/1291-1072/1472) -59.7 E06-W04 N20-N33
      4 1484 1230 (N15W31 L=336)   270     79        0 (1463/1192-1513/1255)  36.9 W29-W33 N12-N16
      5  240 1112 (N10E58 L=246)   489     69        1 ( 226/1075- 258/1147) -68.2 E60-E56 N07-N12
      6 1749  852 (S10W51 L=355)   612     48        1 (1726/ 815-1772/ 879) -53.1 W49-W53 S12-S08
      7  815  817 (S08E14 L=290)   777    183        2 ( 736/ 775- 872/ 861)  26.6 E19-E10 S10-S05
      8 1130  715 (S15W06 L=310)   227     71        0 (1120/ 679-1143/ 753)  90.0 W05-W06 S17-S12
      9  342  650 (S19E53 L=251)   593     84        1 ( 315/ 615- 381/ 676)  26.6 E56-E48 S21-S16
     10 1564  601 (S26W38 L=342)   492     96        1 (1532/ 570-1602/ 639) -26.6 W34-W42 S27-S22
     11 1180  147 (S68W18 L=323)   721     43        1 (1162/ 118-1196/ 176) -43.0 W13-W30 S76-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