Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  833 1741 (N32E33 L=314)   594     85        1 ( 795/1699- 864/1774) -26.6 E37-E29 N29-N34
      2  986 1591 (N25E17 L=331)   373    114        2 ( 969/1533-1002/1645) -54.5 E18-E13 N21-N28
      3  852 1539 (N19E23 L=325)   488    186        1 ( 766/1515- 943/1560)  15.3 E28-E17 N17-N20
      4 1582 1356 (N26W27 L= 14)   598    153        2 (1543/1285-1632/1398)  77.5 W24-W30 N21-N30
      5 1686 1296 (N26W36 L= 24)   412    107        1 (1652/1250-1718/1339) -33.7 W32-W38 N23-N27
      6 1299 1292 (N14W09 L=357)   385    127        3 (1270/1241-1339/1328) -68.2 W07-W12 N11-N16
      7  944 1231 (N02E09 L=339)   580    241        0 ( 827/1208-1049/1264)  26.6 E15-E03 N00-N06
      8  313 1205 (S11E50 L=298)   395    112        2 ( 289/1152- 342/1251)  66.8 E51-E48 S14-S07
      9 1607 1203 (N17W31 L= 19)   425    124        2 (1566/1166-1667/1237)   7.1 W27-W35 N15-N19
     10 1842  994 (N13W55 L= 43)   756    140        0 (1831/ 920-1855/1067) -71.6 W53-W57 N08-N17
     11  206  977 (S26E58 L=290)  1547    185        3 ( 130/ 934- 277/1010)  14.7 E69-E50 S28-S23
     12 1839  805 (N03W60 L= 48)   407     40        2 (1815/ 781-1859/ 829)   0.0 W57-W63 N01-N04
     13 1461  731 (S12W33 L= 21)   227     53        1 (1435/ 707-1490/ 751)  33.7 W31-W34 S13-S10
     14  198  668 (S45E68 L=280)  1236    117        2 ( 137/ 639- 257/ 690)  -6.6 E85-E57 S46-S43
     15 1151  649 (S25W17 L=  5)   377     86        0 (1109/ 629-1199/ 671)   9.5 W14-W20 S26-S23
     16 1465  515 (S24W43 L= 31)  2104    520        3 (1252/ 459-1712/ 599)  38.9 W30-W58 S32-S11
     17 1767  487 (S15W73 L= 61)   832     72        2 (1732/ 471-1806/ 506)  12.1 W67-W81 S16-S13
  

    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