Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  874 1601 (N26E26 L=223)  1843    356        4 ( 759/1484-1021/1682)  63.4 E28-E18 N16-N34
      2 1077 1581 (N30E12 L=236)  1318    178        2 (1014/1507-1128/1647) -31.0 E17-E07 N25-N34
      3 1529 1303 (N25W24 L=273)  2795    565        5 (1391/1060-1596/1480) -52.3 W11-W31 N11-N35
      4  205 1217 (S11E63 L=186)   851     82        3 ( 171/1183- 235/1271)  80.5 E67-E58 S13-S07
      5  160 1094 (S19E66 L=182)  2394    257        2 (  82/1015- 255/1186)  22.5 E85-E55 S24-S12
      6 1800  953 (N13W53 L=302)   510     95        0 (1782/ 907-1813/1006)  90.0 W51-W54 N09-N15
      7  427  955 (S22E36 L=213)   945    276        0 ( 366/ 862- 539/1005)  42.5 E40-E28 S28-S16
      8  884  838 (S18E03 L=245)   392    181        0 ( 850/ 776- 934/ 879) -80.5 E05-W00 S20-S14
      9  457  845 (S28E32 L=216)   683    128        2 ( 396/ 829- 532/ 864)  15.3 E37-E26 S29-S24
     10 1351  768 (S10W25 L=274)   415    149        0 (1287/ 724-1408/ 793)   6.3 W21-W30 S11-S08
     11  766  377 (S51W02 L=251)  1595    263        4 ( 673/ 299- 860/ 467) -15.3 E08-W14 S54-S46
     12  963  238 (S55W28 L=277)   341     45        1 ( 937/ 217- 987/ 253)  33.7 W25-W31 S56-S53
  

    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