Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1426 1687 (N31W50 L=262)   745    168        0 (1349/1664-1506/1706) -31.0 W43-W55 N27-N34
      2  910 1509 (N35W06 L=218)  1677    639        0 ( 642/1451-1197/1555) -16.1 E14-W24 N27-N38
      3 1473 1373 (N12W36 L=248)   709    152        0 (1425/1314-1524/1437) -77.5 W34-W37 N07-N16
      4 1899 1170 (S11W67 L=280)  1197    140        0 (1864/1129-1959/1216) -21.3 W63-W81 S14-S07
      5  404  871 (N08E42 L=170)  1722    522        2 ( 234/ 706- 565/1005) -45.0 E55-E34 S04-N18
      6  216  824 (N09E61 L=151)   512    129        0 ( 186/ 769- 253/ 878) -66.8 E64-E58 N05-N12
      7 1394  749 (S21W15 L=228)   651    164        0 (1351/ 686-1425/ 825)  40.6 W11-W18 S23-S17
      8 1755  649 (S38W46 L=258)  2844    387        3 (1647/ 520-1854/ 810)  32.5 W34-W56 S45-S30
      9  527  408 (S21E51 L=161)  1068    228        1 ( 495/ 319- 563/ 506)  49.8 E56-E43 S28-S15
     10  685  415 (S24E39 L=173)   544    157        0 ( 665/ 339- 714/ 481)  42.0 E43-E33 S29-S20
     11  820  392 (S29E31 L=181)   355     93        0 ( 796/ 350- 845/ 440)  35.5 E34-E27 S31-S26
     12 1270  299 (S50E05 L=208)  1467    287        2 (1153/ 249-1409/ 362)   0.0 E16-W09 S51-S47
  

    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