Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  947 1779 (N37E33 L=145)   301     58        1 ( 929/1750- 963/1808) -51.3 E35-E30 N34-N39
      2 1646 1517 (N44W29 L=207)   687    106        3 (1600/1476-1695/1546)  -5.7 W24-W35 N42-N46
      3  274 1406 (S03E63 L=115)   541     72        0 ( 256/1376- 296/1436)  76.0 E65-E60 S04-N00
      4  259 1317 (S08E60 L=117)   368     76        1 ( 232/1297- 289/1339)  36.9 E64-E57 S09-S06
      5  897 1241 (N02E13 L=165)   397    105        1 ( 840/1225- 946/1261)  38.7 E15-E10 N00-N04
      6 1806 1173 (N27W48 L=226)   343     85        0 (1775/1151-1848/1189)  32.0 W45-W53 N24-N29
      7  205 1075 (S23E60 L=118)   765    148        1 ( 143/1021- 257/1120)  26.6 E69-E54 S27-S19
      8  350  859 (S33E42 L=136)   659    168        2 ( 320/ 781- 395/ 926)  59.0 E44-E38 S38-S28
      9  458  841 (S32E32 L=146)   328     66        1 ( 436/ 809- 494/ 874) -21.8 E33-E28 S33-S29
     10  546  793 (S32E24 L=154)   538     97        1 ( 499/ 771- 593/ 824)  35.5 E27-E20 S34-S29
     11 1055  789 (S18W08 L=186)   214    113        0 (1005/ 777-1100/ 818)  45.0 W05-W10 S20-S15
     12 1128  686 (S22W16 L=194)   780    276        0 (1094/ 598-1202/ 760) -84.8 W13-W21 S26-S15
  

    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