Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1139 1836 (N49E22 L= 96)   502    105        1 (1103/1801-1185/1855)  -5.7 E25-E15 N46-N51
      2  945 1634 (N27E22 L= 95)  1463    364        3 ( 902/1508-1011/1777) -41.6 E31-E13 N18-N37
      3  597 1726 (N25E52 L= 66)   842    159        3 ( 545/1678- 659/1763)  -8.1 E59-E45 N21-N27
      4  395 1599 (N13E61 L= 56)   564    113        1 ( 376/1546- 411/1654) -53.1 E64-E58 N08-N16
      5  526 1563 (N12E47 L= 71)   475    135        0 ( 505/1497- 556/1615) -63.4 E48-E43 N08-N16
      6 1361 1395 (N23W10 L=127)   918    234        1 (1327/1297-1382/1501) -54.0 W05-W13 N16-N29
      7 1688 1441 (N37W35 L=152)   549    124        0 (1641/1411-1736/1491) -11.3 W30-W40 N36-N39
      8  554 1387 (N02E37 L= 80)   698    129        3 ( 487/1358- 607/1417)  14.0 E42-E34 N00-N04
      9  717  733 (S31E11 L=106)   312     98        1 ( 669/ 712- 766/ 755)   8.1 E14-E07 S32-S29
     10  835  585 (S37W01 L=118)  1019    330        0 ( 669/ 566- 974/ 607)  28.8 E11-W09 S43-S32
     11 1117  482 (S35W23 L=141)   987    248        2 (1063/ 413-1196/ 587) -20.6 W15-W31 S38-S30
     12 1577  409 (S24W61 L=178)   431    106        0 (1552/ 387-1625/ 443)  23.2 W56-W65 S25-S22
  

    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