Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1459 1591 (N24W47 L=179)   545    122        2 (1402/1561-1515/1617) -50.2 W43-W49 N20-N27
      2 1311 1450 (N21W30 L=162)  1188    502        2 (1096/1396-1511/1519) -36.5 W19-W42 N14-N31
      3 1783 1311 (S01W59 L=190)   642     95        2 (1749/1271-1816/1341) -45.0 W55-W61 S03-N01
      4 1083 1216 (N15W08 L=140)   319    164        0 (1005/1192-1156/1231) -18.4 W03-W12 N12-N16
      5  747 1167 (N21E12 L=119)   352    115        2 ( 697/1139- 797/1182) -26.6 E15-E09 N18-N22
      6 1880 1024 (S20W60 L=192)  1220    185        2 (1854/ 940-1907/1111)  68.2 W56-W64 S24-S14
      7  280 1078 (N26E50 L= 82)   492     88        3 ( 247/1050- 328/1102) -26.6 E53-E45 N24-N28
      8  178 1055 (N27E63 L= 69)   724     88        2 ( 150/1020- 214/1096) -45.0 E67-E58 N24-N29
      9  484  988 (N17E34 L= 98)   218     90        0 ( 441/ 980- 531/ 997) -15.9 E37-E30 N15-N17
     10 1574  767 (S24W27 L=158)   239     43        1 (1545/ 744-1599/ 782) -14.0 W24-W28 S26-S23
     11 1253  610 (S24W02 L=134)   956    226        5 (1221/ 509-1271/ 731)  65.2 E01-W05 S30-S17
     12  477  497 (S11E51 L= 81)   411     65        1 ( 453/ 467- 506/ 533) -90.0 E53-E49 S13-S08
     13 1211  455 (S33E05 L=126)   300     40        1 (1185/ 432-1237/ 473) -33.7 E06-E03 S34-S31
  

    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