Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  880 1794 (N58W24 L=343)  3902    466        7 ( 702/1660-1045/1926)   7.4 E01-W62 N53-N62
      2 1454 1724 (N31W58 L= 18)   418     42        0 (1432/1699-1476/1751)  53.1 W55-W61 N29-N33
      3 1482 1613 (N24W50 L= 10)   386    118        0 (1438/1597-1527/1637) -14.0 W47-W55 N23-N26
      4  381 1535 (N56E42 L=277)   397     47        0 ( 359/1509- 405/1557)  45.0 E46-E38 N53-N57
      5  966 1441 (N31W09 L=328)   412    120        0 ( 918/1402- 992/1491)  33.7 W05-W11 N29-N33
      6  266 1264 (N39E52 L=268)   523     69        2 ( 224/1242- 302/1277) -11.3 E57-E47 N36-N39
      7 1622  879 (S19W32 L=352)   559    107        0 (1605/ 823-1637/ 943)  60.3 W30-W34 S22-S15
      8 1461  755 (S21W18 L=338)   344    130        0 (1429/ 700-1490/ 814)  39.8 W15-W21 S23-S18
      9 1530  503 (S41W19 L=338)  1591    271        2 (1474/ 439-1621/ 584) -12.1 W13-W27 S45-S34
     10 1107  490 (S27E11 L=309)   342    110        0 (1060/ 460-1160/ 524) -60.3 E12-E08 S30-S23
     11  939  460 (S24E22 L=298)   265    101        0 ( 892/ 440- 995/ 486)   0.0 E25-E17 S25-S23
     12  746  315 (S29E42 L=277)  1272    360        1 ( 604/ 266- 880/ 391) -42.9 E48-E34 S35-S21
  

    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