Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  568 1612 (N15E47 L=280)   409    122        0 ( 553/1566- 590/1667) -60.3 E48-E44 N12-N19
      2 1305 1532 (N31W02 L=329)   414     99        1 (1262/1504-1356/1554)   0.0 E01-W05 N29-N32
      3 1707 1505 (N45W35 L=  2)   759    133        1 (1683/1454-1724/1550) -90.0 W33-W37 N40-N49
      4  517 1405 (N02E41 L=286)   399     90        2 ( 498/1366- 541/1448) -45.0 E43-E38 N00-N04
      5 1543 1410 (N31W21 L=348)   262     47        0 (1517/1385-1566/1433) -11.3 W18-W23 N30-N32
      6  702 1381 (N05E28 L=299)   233     95        1 ( 663/1362- 745/1395)  36.9 E30-E25 N02-N05
      7  276 1205 (S14E54 L=273)  1786    407        3 ( 160/1175- 385/1225)   5.0 E66-E43 S16-S11
      8 1831  982 (N16W53 L= 20)   532     99        0 (1809/ 931-1843/1021) -45.0 W50-W55 N13-N18
      9  838  856 (S21E06 L=321)   259    109        0 ( 823/ 807- 852/ 918) -50.2 E08-E03 S23-S17
     10 1769  831 (N06W53 L= 20)   795    184        1 (1725/ 793-1820/ 890)  -9.5 W46-W59 N04-N07
     11 1361  781 (S10W26 L=353)   211     89        0 (1345/ 734-1377/ 825) -68.2 W25-W27 S12-S07
     12 1412  651 (S15W34 L=  1)   348    160        0 (1349/ 601-1454/ 704) -11.3 W28-W38 S17-S14
     13 1520  681 (S10W40 L=  7)   335     86        3 (1486/ 658-1552/ 703)  36.9 W37-W41 S12-S08
     14 1697  616 (S08W56 L= 23)   359     47        0 (1676/ 584-1720/ 641)  90.0 W54-W57 S09-S05
     15 1597  587 (S13W49 L= 16)   407    108        0 (1548/ 564-1635/ 630)  56.3 W46-W51 S15-S09
     16  996  561 (S33W13 L=340)  1196    271        2 ( 860/ 524-1116/ 625)  42.7 W04-W18 S38-S26
     17 1444  407 (S27W49 L= 16)  1843    361        3 (1323/ 328-1545/ 481)   2.3 W36-W61 S28-S24
  

    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