Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  841 1862 (N55E47 L=208)   934    138        0 ( 781/1825- 900/1884)  -4.8 E59-E35 N52-N57
      2 1150 1735 (N50E06 L=249)   347     63        1 (1129/1704-1179/1754)  53.1 E07-E02 N47-N52
      3  468 1680 (N32E63 L=192)   871     87        3 ( 443/1645- 501/1719) -13.2 E73-E56 N29-N34
      4 1300 1634 (N45W10 L=265)   405     43        0 (1276/1606-1325/1656) -11.3 W07-W12 N42-N46
      5 1795 1357 (N34W59 L=313)   662     85        0 (1755/1341-1832/1378)   9.5 W52-W64 N31-N35
      6 1587 1241 (N22W35 L=289)   518    172        1 (1526/1184-1643/1297) -26.6 W29-W39 N19-N24
      7  282  940 (S17E51 L=203)   387     66        0 ( 257/ 921- 318/ 963)  -8.1 E54-E47 S18-S15
      8  266  877 (S21E53 L=202)   447     87        2 ( 235/ 853- 287/ 912)  26.6 E57-E50 S23-S18
      9 1453  441 (S29W44 L=298)   986    175        1 (1370/ 422-1536/ 471)  24.8 W37-W50 S32-S26
     10 1472  378 (S33W50 L=305)   433     61        3 (1436/ 357-1501/ 402)  36.9 W47-W52 S35-S31
     11 1043  196 (S59W28 L=283)   537     58        0 (1019/ 169-1063/ 227) -45.0 W24-W32 S62-S55
  

    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