Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  866 1793 (N40E32 L=329)   844    127        3 ( 805/1761- 915/1819)   0.0 E38-E27 N37-N42
      2 1103 1716 (N39E11 L=350)   578    141        1 (1045/1671-1145/1748) -12.5 E15-E06 N36-N41
      3 1216 1696 (N40E02 L=359)   968    125        2 (1174/1644-1257/1748) -36.9 E06-W02 N36-N43
      4  494 1668 (N24E54 L=307)   909    184        5 ( 445/1602- 550/1717)  90.0 E57-E49 N19-N28
      5  384 1619 (N20E62 L=299)  1115    206        4 ( 326/1536- 448/1685)  80.5 E65-E57 N13-N25
      6 1161 1358 (N15W02 L=  2)   211     56        0 (1131/1337-1193/1374) -18.4 W00-W03 N14-N16
      7 1567 1329 (N22W29 L= 29)   256     62        2 (1541/1303-1592/1359) -14.0 W26-W30 N21-N23
      8 1096 1268 (N09E00 L=  0)   186     78        0 (1061/1239-1133/1287)  56.3 E02-W01 N06-N10
      9  245 1161 (S10E57 L=304)   580     84        4 ( 219/1118- 272/1203)  78.7 E59-E54 S12-S07
     10  174 1110 (S14E64 L=296)   442     97        0 ( 155/1071- 204/1146)  39.8 E66-E60 S16-S11
     11  334  989 (S20E46 L=314)   571    212        3 ( 310/ 892- 364/1083)  85.2 E47-E43 S25-S13
     12  280  758 (S34E53 L=308)   494     97        3 ( 241/ 726- 313/ 798)  29.1 E58-E48 S36-S31
     13 1706  671 (S12W54 L= 54)   429    110        1 (1683/ 631-1726/ 728) -39.8 W49-W56 S14-S09
     14 1327  611 (S24W28 L= 29)   751    214        3 (1289/ 511-1371/ 713) -72.9 W25-W31 S30-S17
     15 1555  475 (S26W52 L= 52)  1573    331        3 (1450/ 397-1685/ 580)  63.4 W46-W57 S33-S17
     16 1605  386 (S30W65 L= 65)  1128    168        2 (1532/ 342-1671/ 431)  56.3 W60-W69 S34-S25
  

    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