Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1440 1473 (N18W39 L=327)   270    100        1 (1403/1446-1469/1502)   9.5 W35-W41 N17-N19
      2  610 1387 (N37E17 L=272)   255     62        1 ( 580/1371- 647/1403)  -9.5 E19-E13 N36-N38
      3 1224 1346 (N18W20 L=309)   200     71        0 (1195/1331-1264/1371) -45.0 W18-W22 N15-N19
      4 1200 1162 (N08W13 L=302)   544    137        3 (1122/1141-1268/1182) -24.0 W08-W17 N05-N10
      5 1320 1139 (N04W19 L=308)   377    118        1 (1280/1081-1373/1173)  15.9 W16-W23 N01-N05
      6  825 1071 (N13E10 L=279)   407    167        1 ( 741/1020- 900/1112)   0.0 E15-E04 N11-N14
      7 1797  954 (S21W48 L=336)   343     75        0 (1782/ 909-1812/ 982)  45.0 W45-W49 S23-S19
      8  193  879 (N16E61 L=227)  2454    462        4 (  89/ 835- 306/ 933) -16.9 E74-E50 N11-N20
      9 1363  661 (S25W10 L=299)  1505    337        3 (1310/ 516-1429/ 802)  47.0 W02-W17 S32-S17
     10  752  516 (S16E30 L=258)   293     38        0 ( 725/ 500- 782/ 535) -56.3 E31-E28 S17-S14
     11  653  419 (S20E41 L=247)   328    123        0 ( 603/ 401- 723/ 432) -29.7 E44-E37 S22-S18
     12  595  353 (S23E49 L=239)   342     78        1 ( 566/ 327- 622/ 376) -68.2 E51-E46 S25-S20
  

    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