Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1659 1366 (N25W38 L=197)   954    107        1 (1606/1320-1697/1404) -23.2 W33-W41 N22-N28
      2  909 1291 (N08E11 L=148)   264     87        0 ( 875/1260- 946/1331) -45.0 E13-E07 N06-N10
      3 1724 1217 (N17W44 L=203)   740    103        1 (1682/1183-1771/1258)  -7.1 W40-W48 N15-N19
      4 1131 1128 (N01W05 L=163)   382    116        2 (1085/1088-1184/1170) -26.6 W01-W08 N00-N03
      5 1550 1072 (N05W32 L=190)   229     61        0 (1524/1048-1579/1099) -14.0 W29-W33 N03-N06
      6 1151 1009 (S05W08 L=166)   326    165        0 (1071/ 977-1224/1032)   0.0 W02-W12 S06-S04
      7 1213  891 (S11W13 L=172)   399    188        0 (1182/ 811-1266/ 957)  78.7 W11-W15 S16-S06
      8  227  882 (S23E58 L=101)   647     88        1 ( 203/ 851- 265/ 924) -76.0 E60-E53 S24-S20
      9  168  871 (S24E66 L= 93)   686     83        2 ( 151/ 837- 189/ 917) -63.4 E68-E63 S25-S21
     10 1849  748 (S07W66 L=225)   556     43        2 (1823/ 724-1870/ 776)  36.9 W62-W69 S08-S05
     11 1243  683 (S23W19 L=178)   344     83        1 (1193/ 668-1285/ 701)   0.0 W15-W22 S23-S21
     12 1749  612 (S17W61 L=220)   789    122        1 (1690/ 584-1806/ 637)  20.0 W56-W67 S20-S14
  

    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