Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  650 1681 (N46E37 L=129)  2996    464        4 ( 528/1603- 801/1734)   4.8 E46-E21 N39-N50
      2 1271 1636 (N45W19 L=184)   405     78        1 (1252/1595-1289/1669)  80.5 W17-W21 N41-N47
      3  983 1617 (N42E06 L=160)   366     72        0 ( 959/1589-1007/1646)  45.0 E07-E03 N39-N44
      4 1399 1539 (N37W28 L=193)   294     59        0 (1372/1522-1431/1558)  11.3 W25-W30 N36-N38
      5  961 1409 (N27E06 L=160)   261     80        1 ( 921/1382-1000/1433)  -9.5 E08-E02 N25-N28
      6  508 1233 (N14E35 L=130)   573    119        1 ( 462/1192- 537/1262) -23.2 E39-E32 N11-N15
      7  583 1249 (N15E30 L=135)   250     91        0 ( 542/1235- 625/1260)   0.0 E32-E26 N14-N15
      8  621 1192 (N12E27 L=139)   315     48        1 ( 604/1158- 643/1230) -78.7 E28-E24 N09-N14
      9  498 1152 (N09E35 L=130)   704    234        2 ( 433/1110- 610/1204) -15.9 E40-E26 N06-N11
     10 1587  575 (S25W43 L=208)   502    124        2 (1548/ 543-1627/ 605) -26.6 W38-W46 S27-S22
     11 1406  419 (S37W32 L=198)  1183    214        1 (1341/ 393-1490/ 469)   9.5 W26-W39 S38-S32
     12  986  389 (S40E01 L=164)  1085    264        2 ( 857/ 352-1127/ 421)  13.4 E11-W10 S43-S36
  

    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