Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  377 1443 (N29E44 L=171)   471     86        1 ( 350/1400- 399/1482) -45.0 E48-E41 N26-N32
      2 1470 1373 (N16W31 L=246)   334    116        1 (1444/1323-1516/1409)  38.7 W28-W35 N13-N18
      3  123 1148 (N13E69 L=147)   876    134        3 (  95/1107- 160/1191)  33.7 E74-E62 N10-N15
      4 1108  939 (S07W04 L=220)   224    124        0 (1094/ 876-1124/ 996)  74.1 W03-W05 S11-S04
      5  904  952 (S05E08 L=208)   440     63        1 ( 865/ 918- 931/ 984) -76.0 E09-E05 S07-S03
      6  197  883 (S03E60 L=155)   349     58        0 ( 186/ 854- 214/ 916)  45.0 E61-E58 S05-S01
      7  312  700 (S15E53 L=163)  4983    612       13 ( 195/ 545- 442/ 844)  40.4 E68-E39 S24-S07
      8 1404  548 (S34W23 L=239)   553    166        0 (1317/ 511-1473/ 587)  18.4 W16-W28 S35-S31
      9  946  518 (S32E10 L=206)   585    125        2 ( 889/ 502-1024/ 539) -12.5 E13-E04 S33-S30
     10 1139  449 (S39W03 L=219)   602    163        1 (1067/ 426-1228/ 472)  -5.2 E01-W10 S40-S36
     11  416  439 (S33E56 L=160)   830     89        2 ( 373/ 412- 460/ 470)   4.1 E63-E49 S34-S30
  

    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