Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1236 1679 (N45E02 L=309)   304    106        1 (1188/1656-1281/1704)  -7.1 E06-W02 N43-N46
      2  454 1579 (N19E55 L=256)   570     84        3 ( 422/1539- 483/1619) -80.5 E56-E53 N16-N22
      3  707 1554 (N23E34 L=278)   284      0       -1 ( 689/1530- 728/1579) -14.0 E35-E31 N21-N24
      4 1038 1507 (N27E11 L=301)   934    345        1 ( 897/1412-1160/1562)  -8.1 E21-W00 N24-N29
      5 1249 1458 (N29W04 L=316)   297    135        2 (1208/1405-1284/1509) -35.5 W00-W07 N26-N31
      6  567 1374 (N09E37 L=275)   648    152        3 ( 510/1317- 612/1428)  63.4 E38-E34 N04-N12
      7 1784 1213 (N27W51 L=  3)  1518    217        4 (1752/1111-1807/1319)  81.9 W47-W55 N19-N34
      8  311 1109 (S11E50 L=262)   674    101        3 ( 278/1060- 348/1138)  63.4 E52-E46 S14-S08
      9 1627 1016 (N11W38 L=350)   402    100        0 (1599/ 966-1657/1052)  63.4 W36-W40 N07-N13
     10 1888  942 (N13W67 L= 19)  1266    139        3 (1862/ 874-1907/ 997)  81.9 W64-W71 N08-N15
     11  364  920 (S22E42 L=269)   807    113        2 ( 324/ 884- 404/ 952) -26.6 E46-E38 S23-S19
     12 1825  912 (N09W59 L= 10)   468     74        2 (1805/ 880-1842/ 950)  78.7 W57-W62 N06-N11
     13 1289  605 (S21W27 L=338)  1240    328        4 (1233/ 480-1384/ 695) -86.2 W24-W31 S28-S13
     14  359  590 (S44E46 L=266)   556     69        1 ( 341/ 553- 385/ 622)  78.7 E48-E42 S47-S41
     15  786  583 (S34E05 L=306)   378    148        1 ( 722/ 561- 855/ 608)  24.0 E09-W00 S36-S31
     16 1396  479 (S26W39 L=351)   434    195        1 (1299/ 459-1480/ 502)   7.6 W31-W46 S27-S24
  

    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