Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1471 1335 (N21W32 L=347)  2374    644        9 (1200/1225-1642/1436) -30.4 W14-W43 N12-N29
      2 1827 1281 (N14W64 L= 18)  1022    119        3 (1794/1232-1859/1327) -71.6 W59-W67 N10-N16
      3  261 1108 (N11E55 L=260)  1275    263        4 ( 210/1066- 312/1180)  56.3 E60-E49 N08-N15
      4  521 1054 (N08E33 L=282)   256     48        0 ( 493/1037- 539/1075)  33.7 E34-E31 N07-N09
      5 1860  935 (S07W63 L= 17)  1155    249        3 (1800/ 891-1907/ 975) -23.2 W55-W70 S10-S04
      6  627  912 (S01E26 L=289)   233     55        1 ( 600/ 888- 655/ 931)  14.0 E27-E23 S02-N00
      7 1512  750 (S16W31 L=345)  1167    295        1 (1402/ 638-1610/ 891)  45.0 W23-W38 S22-S07
      8 1606  553 (S30W41 L=356)  1658    300        4 (1502/ 513-1701/ 598)  10.0 W32-W49 S32-S27
      9  408  551 (S24E49 L=266)   418     46        0 ( 387/ 522- 426/ 576) -76.0 E49-E46 S26-S22
  

    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