Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1516 1628 (N31W44 L=203)   387    125        1 (1493/1582-1544/1658)  33.7 W41-W47 N28-N33
      2 1234 1605 (N33W23 L=181)   274     61        1 (1212/1572-1251/1634)  63.4 W21-W24 N30-N35
      3 1720 1417 (N15W54 L=212)   476    130        0 (1691/1376-1762/1442)  29.7 W50-W58 N12-N16
      4  227 1422 (N34E57 L=101)   645    136        0 ( 196/1405- 267/1436) -12.5 E61-E52 N32-N34
      5  177 1143 (N16E60 L= 98)  3430    553       10 (  86/1023- 242/1276)  33.0 E74-E52 N10-N23
      6  446 1146 (N13E35 L=123)   870    282        1 ( 369/1041- 515/1243) -55.3 E41-E31 N06-N19
      7  649  995 (N02E23 L=135)   546    114        2 ( 617/ 931- 675/1056)  56.3 E25-E20 S01-N05
      8 1724  956 (S12W45 L=204)   370    127        1 (1674/ 927-1764/ 995)  23.2 W41-W49 S13-S10
      9 1074  642 (S24E01 L=157)   462    130        0 (1006/ 625-1151/ 666)   0.0 E05-W03 S25-S22
     10  342  553 (S21E56 L=103)   742    178        1 ( 322/ 495- 361/ 614)  41.2 E59-E51 S24-S17
     11 1553  526 (S38W34 L=192)   320    115        1 (1531/ 483-1571/ 573) -60.3 W31-W36 S41-S34
     12  972  395 (S40E13 L=145)   796    161        1 ( 908/ 351-1036/ 452)  22.6 E19-E07 S42-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