Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  837 1763 (N44E21 L= 62)   537     77        1 ( 796/1742- 879/1784)  -8.1 E24-E17 N42-N45
      2  923 1712 (N40E13 L= 71)   318     43        0 ( 907/1686- 943/1739)  76.0 E13-E11 N37-N42
      3 1137 1642 (N36W04 L= 87)  2206    432        2 ( 954/1552-1345/1722) -17.2 E10-W19 N30-N41
      4  556 1627 (N32E39 L= 45)  1102    298        0 ( 475/1549- 617/1696) -32.9 E48-E31 N26-N37
      5 1495 1518 (N30W30 L=113)   590     85        1 (1450/1497-1538/1535)   0.0 W26-W33 N27-N31
      6 1663 1368 (N21W42 L=125)   310     73        0 (1636/1327-1689/1397) -45.0 W39-W43 N18-N23
      7 1040 1337 (N14E01 L= 82)   368     95        0 (1024/1291-1057/1393) -80.5 E01-W00 N11-N17
      8 1265 1313 (N14W13 L= 96)   275     84        0 (1218/1278-1297/1335) -31.0 W10-W15 N12-N15
      9 1643  588 (S26W48 L=131)   320     76        0 (1607/ 573-1680/ 602)  21.8 W45-W50 S27-S24
     10 1441  507 (S34W34 L=117)   293     39        0 (1415/ 490-1471/ 525)   0.0 W31-W36 S34-S32
  

    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