Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  742 1714 (N31E36 L= 65)  1453    265        2 ( 635/1644- 836/1795) -17.7 E51-E26 N26-N35
      2  428 1652 (N21E60 L= 42)  2402    238        4 ( 386/1540- 490/1728) -74.1 E68-E52 N13-N27
      3 1124 1651 (N34E08 L= 94)   704     88        2 (1085/1605-1157/1687) -23.2 E11-E04 N30-N36
      4  731 1590 (N22E31 L= 70)   270     96        0 ( 695/1559- 774/1614)  53.1 E32-E29 N18-N23
      5 1095 1299 (N10E01 L=100)  1304    447        3 (1056/1107-1131/1490) -62.4 E07-W04 N00-N21
      6 1351 1370 (N20W13 L=114)   775    266        2 (1251/1295-1439/1432) -12.1 W05-W19 N17-N22
      7 1772 1292 (N27W47 L=148)   962    156        0 (1744/1234-1799/1366) -76.0 W43-W49 N23-N31
      8 1828 1346 (N32W54 L=156)   482     83        0 (1791/1313-1857/1364)   0.0 W49-W58 N31-N33
      9 1329 1287 (N14W13 L=114)   341     87        3 (1305/1249-1359/1321) -53.1 W10-W15 N12-N16
     10 1820 1116 (N17W52 L=154)   427     47        1 (1797/1091-1840/1144)  71.6 W49-W54 N15-N18
     11 1397  801 (S12W27 L=129)  1067    344        3 (1278/ 750-1543/ 856)   9.5 W18-W36 S14-S07
     12 1212  781 (S17W16 L=118)   354    182        0 (1190/ 733-1269/ 838) -45.0 W14-W20 S19-S13
     13  957  727 (S25W02 L=104)   494    146        2 ( 896/ 677-1011/ 762) -14.0 E02-W06 S27-S23
     14 1714  651 (S12W56 L=157)   410     57        1 (1690/ 631-1738/ 673)   0.0 W52-W58 S13-S10
     15  883  601 (S35W00 L=102)   234     49        1 ( 851/ 588- 917/ 618)   0.0 E02-W02 S35-S33
     16 1479  422 (S31W50 L=151)   597    128        2 (1419/ 399-1535/ 451)   9.5 W43-W55 S32-S29
     17  838  201 (S68W14 L=116)   957    136        3 ( 771/ 166- 901/ 238)   0.0 W03-W27 S70-S65
  

    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