Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  955 1612 (N39E12 L=126)  1073    219        3 ( 876/1545-1041/1659) -21.8 E19-E04 N34-N41
      2  474 1513 (N27E47 L= 91)  1629    245        4 ( 409/1422- 532/1589) -32.0 E57-E39 N20-N31
      3 1022 1383 (N24E04 L=135)   460    159        2 (1003/1314-1042/1465)  90.0 E05-E02 N19-N29
      4  706 1387 (N21E25 L=114)   258     73        0 ( 669/1365- 745/1407)  26.6 E27-E22 N19-N22
      5  616  777 (S17E24 L=114)  2769    615        5 ( 399/ 588- 797/ 991)  43.1 E44-E13 S32-S02
      6 1861  728 (S10W71 L=209)   879     93        3 (1829/ 688-1891/ 769) -68.2 W66-W76 S12-S07
      7 1647  683 (S14W47 L=185)   830    145        2 (1589/ 659-1723/ 706)   0.0 W41-W54 S15-S12
      8 1570  505 (S27W47 L=186)  1230    252        1 (1476/ 433-1651/ 598)  50.7 W41-W51 S32-S19
      9  286  515 (S40E66 L= 72)   711     68        0 ( 268/ 481- 314/ 546)  18.4 E77-E58 S43-S37
     10 1384  429 (S34W34 L=173)   335     64        0 (1350/ 415-1417/ 443)  -9.5 W31-W37 S34-S32
  

    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