Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  628 1717 (N61E10 L=126)   422     81        0 ( 599/1686- 665/1747)  18.4 E13-E04 N59-N63
      2  552 1507 (N46E21 L=115)   340     55        1 ( 531/1478- 579/1529)  21.8 E22-E17 N44-N47
      3  996 1486 (N32W12 L=147)   292     55        1 ( 967/1465-1025/1510) -76.0 W09-W13 N29-N33
      4 1207 1423 (N22W22 L=158)   206     97        0 (1166/1397-1251/1439)   0.0 W18-W25 N21-N22
      5  460 1154 (N26E32 L=103)   539    149        0 ( 442/1087- 487/1220)  63.4 E34-E29 N21-N29
      6 1819  975 (S20W52 L=188)   699    157        0 (1780/ 909-1872/1022)  14.0 W46-W59 S22-S18
      7 1750  748 (S32W43 L=179)   727    201        0 (1674/ 697-1823/ 819)  10.0 W34-W51 S33-S29
      8  939  675 (S13E14 L=122)   204     41        1 ( 922/ 648- 957/ 709)  63.4 E15-E12 S14-S10
      9  573  569 (S10E40 L= 96)   455     88        0 ( 555/ 519- 599/ 605)  39.8 E42-E36 S13-S08
     10 1461  512 (S38W15 L=150)   722    258        2 (1355/ 457-1577/ 574)   0.0 W05-W25 S39-S36
     11  648  342 (S26E47 L= 89)  1477    178        2 ( 600/ 293- 725/ 421) -77.5 E52-E41 S29-S20
     12  531  319 (S25E59 L= 77)   447     79        0 ( 493/ 308- 571/ 336)   5.2 E64-E53 S26-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