Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1495 1723 (N30W59 L=128)  1614    242        3 (1425/1622-1564/1798)  69.0 W53-W65 N23-N36
      2  678 1759 (N62E02 L= 66)   520     85        1 ( 650/1717- 708/1783)  45.0 E05-W02 N57-N65
      3  381 1641 (N62E49 L= 19)   810     87        1 ( 341/1615- 419/1667)   8.7 E60-E41 N60-N64
      4 1308 1518 (N24W32 L=100)   307     76        0 (1266/1501-1349/1535) -18.4 W28-W34 N23-N26
      5 1104 1423 (N25W16 L= 84)  1152    238        2 (1008/1355-1209/1451) -32.5 W10-W21 N21-N28
      6 1396 1197 (N05W25 L= 94)   327    130        0 (1342/1176-1455/1222) -29.7 W22-W29 N02-N07
      7 1639 1137 (S05W39 L=108)   257     72        1 (1626/1104-1654/1173)  45.0 W37-W41 S06-S03
      8 1869 1042 (S17W58 L=127)   733    127        0 (1848/ 990-1907/1100)  24.0 W54-W64 S20-S15
      9  691 1049 (N14E18 L= 50)   262     79        0 ( 648/1030- 726/1076)   0.0 E21-E15 N13-N14
     10  950  988 (N05E05 L= 63)   258     15        0 ( 926/ 964- 973/1011)   0.0 E06-E03 N03-N05
     11  725  737 (S04E24 L= 44)  1159    371        1 ( 574/ 640- 834/ 855)  15.3 E36-E14 S06-N00
     12  230  727 (N05E62 L=  6)   466     62        0 ( 203/ 701- 253/ 754) -45.0 E64-E59 N03-N07
     13 1510  478 (S42W18 L= 86)  1069    260        1 (1407/ 403-1606/ 551)  14.7 W08-W27 S44-S39
     14  399  459 (S13E59 L=  9)  1125    244        1 ( 355/ 381- 453/ 548)  74.7 E64-E54 S19-S07
     15  475  442 (S16E53 L= 15)   529    108        1 ( 443/ 402- 512/ 489) -90.0 E54-E50 S18-S12
     16 1263  342 (S44E05 L= 63)   375     84        0 (1217/ 326-1313/ 356) -14.0 E08-W00 S45-S41
     17 1000  250 (S42E28 L= 40)   914    175        2 ( 907/ 234-1085/ 273) -12.1 E35-E21 S45-S38
  

    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