Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1388 1389 (N16W32 L= 55)   232     36        0 (1372/1362-1410/1407)  18.4 W30-W33 N15-N17
      2  714 1046 (N16E17 L=  5)  2304    491        7 ( 576/ 914- 891/1238) -65.2 E22-E09 N04-N30
      3 1065 1193 (N15W07 L= 29)   743    272        3 ( 944/1134-1165/1238) -45.0 W00-W11 N09-N20
      4  161 1131 (N32E67 L=315)   559    104        0 ( 135/1083- 189/1171) -26.6 E73-E62 N29-N34
      5  385  994 (N20E42 L=341)   308    142        1 ( 329/ 975- 446/1026) -16.7 E46-E36 N18-N23
      6 1774  836 (S26W45 L= 67)   483     46        1 (1742/ 814-1803/ 863)  -8.1 W41-W48 S27-S24
      7  389  764 (N07E47 L=336)   430     84        0 ( 355/ 737- 420/ 809)  15.9 E50-E43 N05-N08
      8 1541  579 (S35W21 L= 43)  1773    313        4 (1409/ 496-1633/ 682)  11.3 W09-W29 S38-S29
      9  580  370 (S21E50 L=332)   402     72        1 ( 567/ 341- 599/ 406)  45.0 E53-E47 S23-S18
     10 1263  324 (S44E07 L= 15)   657     79        0 (1231/ 288-1301/ 361)  26.6 E10-E02 S46-S42
     11 1193  261 (S47E16 L=  7)   444     39        1 (1173/ 240-1220/ 291)  71.6 E18-E13 S49-S43
     12 1131  228 (S48E23 L=360)   461     64        1 (1105/ 207-1163/ 249) -36.9 E24-E19 S50-S45
  

    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