Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  666 1600 (N17E39 L=358)   262     59        0 ( 637/1574- 698/1618)  63.4 E40-E37 N14-N18
      2 1352 1352 (N21W10 L= 47)   281    110        0 (1307/1300-1398/1403) -29.7 W05-W13 N19-N23
      3 1168 1333 (N15E01 L= 37)   254     61        0 (1134/1313-1197/1361)  56.3 E02-W00 N13-N17
      4 1752 1093 (N20W43 L= 81)  1146    207        0 (1729/ 997-1774/1206) -57.5 W39-W47 N14-N26
      5  911 1137 (S03E09 L= 28)   197     85        0 ( 897/1094- 922/1179) -63.4 E10-E08 S05-S01
      6 1205  964 (S05W11 L= 49)   378     67        1 (1165/ 937-1232/ 995)  63.4 W09-W13 S06-S02
      7 1878  919 (N15W61 L= 98)  1358    142        2 (1850/ 856-1899/ 975) -60.3 W56-W63 N10-N18
      8  534  818 (S31E25 L= 13)  1342    449        2 ( 356/ 766- 720/ 885)   9.8 E40-E11 S33-S26
      9 1365  760 (S11W27 L= 64)   558    152        0 (1327/ 685-1421/ 827)  90.0 W25-W28 S15-S05
     10 1250  611 (S22W26 L= 64)   490     87        1 (1195/ 583-1296/ 640)  51.3 W23-W28 S24-S19
     11 1621  350 (S23W67 L=105)   768     84        0 (1585/ 322-1645/ 384)  90.0 W65-W69 S25-S20
  

    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