Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  633 1526 (N35E23 L= 63)   526     91        0 ( 594/1496- 674/1552) -29.7 E26-E19 N32-N37
      2 1203 1426 (N22W16 L=101)   653    226        0 (1091/1405-1310/1442) -15.9 W08-W22 N19-N23
      3  170 1238 (N21E65 L= 21)  1247    106        2 ( 149/1183- 203/1293) -45.0 E70-E59 N17-N24
      4 1802 1221 (N03W57 L=142)   940     92        4 (1762/1173-1836/1271) -80.5 W52-W60 N00-N06
      5 1593  832 (S18W35 L=121)   684    161        1 (1546/ 754-1624/ 921)  63.4 W31-W37 S22-S12
      6  860  803 (S12E12 L= 73)   228     83        1 ( 815/ 773- 893/ 835)  18.4 E15-E09 S13-S10
      7 1081  678 (S22E00 L= 86)   273    118        0 (1030/ 656-1134/ 712) -38.7 E02-W03 S24-S19
      8  760  585 (S25E22 L= 63)   401     61        1 ( 723/ 569- 807/ 611) -21.8 E24-E19 S26-S23
      9 1093  393 (S42E03 L= 83)  1212    293        1 ( 949/ 334-1225/ 443)   9.9 E15-W08 S45-S39
  

    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