Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  570 1310 (N23E30 L= 15)   294     70        1 ( 548/1279- 599/1342) -45.0 E31-E27 N20-N24
      2 1792 1281 (N16W59 L=104)  1147    146        3 (1735/1228-1828/1336) -54.5 W53-W63 N12-N19
      3  349 1282 (N21E49 L=356)   987    165        1 ( 266/1254- 418/1314)   0.0 E58-E42 N18-N22
      4 1648 1021 (N01W41 L= 86)   386    122        0 (1635/ 957-1662/1083) -90.0 W40-W42 S02-N04
      5  435  969 (N02E39 L=  6)   589    195        1 ( 358/ 915- 507/1018) -24.4 E45-E33 S01-N04
      6 1398  946 (S02W23 L= 68)   235     73        0 (1385/ 905-1416/ 986) -76.0 W22-W24 S04-N00
      7 1819  931 (S06W57 L=102)   768    102        1 (1764/ 914-1872/ 960)  -4.4 W51-W64 S07-S04
      8  963  911 (S03E04 L= 41)   298     72        0 ( 935/ 879-1000/ 945)  45.0 E05-E01 S04-N00
      9  611  846 (S06E27 L= 18)   522    182        0 ( 540/ 814- 656/ 904)  35.5 E31-E23 S07-S02
     10 1297  778 (S12W16 L= 61)  1090    424        1 (1130/ 676-1428/ 897)  32.3 W05-W25 S17-S05
     11  302  711 (S15E54 L=351)   840    211        1 ( 252/ 620- 343/ 798) -65.6 E57-E50 S20-S09
     12  391  655 (S18E47 L=358)   380    119        1 ( 365/ 608- 433/ 695)  33.7 E50-E41 S21-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