Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  871 1867 (N42E47 L=149)  1103    140        2 ( 817/1824- 934/1908)  -9.0 E57-E38 N39-N45
      2 1121 1643 (N32E14 L=182)  1156    223        1 (1059/1551-1173/1728) -32.7 E21-E07 N27-N36
      3  560 1689 (N20E53 L=142)   558     73        2 ( 530/1663- 591/1713)  -7.1 E57-E49 N18-N21
      4  777 1562 (N17E31 L=165)   525    162        1 ( 735/1522- 845/1616) -24.0 E34-E25 N15-N20
      5  710 1501 (N12E32 L=164)   483    134        0 ( 668/1455- 764/1551) -18.4 E36-E27 N10-N14
      6 1242 1382 (N19W02 L=198)   265     95        0 (1217/1340-1273/1429) -38.7 W00-W05 N17-N21
      7  401 1329 (S05E47 L=149)   841    226        0 ( 365/1232- 431/1422) -79.7 E50-E44 S11-N00
      8 1550 1303 (N24W24 L=220)   828    225        0 (1515/1218-1572/1411) -67.4 W19-W26 N18-N30
      9 1709 1310 (N31W37 L=233)   810    193        1 (1639/1245-1748/1382) -31.0 W29-W41 N27-N33
     10  506 1269 (S06E37 L=158)   724    212        0 ( 450/1184- 555/1358)  80.5 E39-E35 S12-N00
     11  982 1286 (N06E09 L=186)   565    150        2 ( 933/1230-1034/1344)  76.0 E10-E07 N02-N10
     12 1561  918 (N03W34 L=230)   313     68        3 (1532/ 893-1595/ 943)  11.3 W31-W36 N01-N04
     13  174  863 (S37E65 L=130)   806    124        0 ( 136/ 826- 220/ 908)  18.4 E75-E57 S39-S33
     14 1352  774 (S11W26 L=222)   205     23        0 (1329/ 758-1371/ 792)  18.4 W24-W27 S11-S09
     15 1264  704 (S17W23 L=219)   219     70        1 (1227/ 676-1289/ 731)  76.0 W22-W24 S19-S15
     16  988  678 (S27W08 L=204)   403     96        2 ( 943/ 654-1035/ 703)   9.5 W05-W11 S28-S25
     17 1465  355 (S29W57 L=253)   930    114        0 (1419/ 333-1530/ 387)  15.3 W51-W63 S30-S26
  

    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