Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1306 1760 (N53W39 L= 55)   830    105        1 (1264/1728-1356/1780)   0.0 W32-W45 N50-N54
      2 1646 1576 (N34W61 L= 77)   512     97        0 (1592/1559-1690/1593) -15.3 W54-W66 N32-N36
      3  547 1485 (N37E35 L=341)   553    122        0 ( 515/1433- 587/1533)  60.3 E37-E32 N32-N39
      4  697 1432 (N32E21 L=355)   933    188        1 ( 640/1378- 778/1501) -39.3 E26-E15 N28-N37
      5  832 1197 (N16E11 L=  5)   747    345        1 ( 709/1099- 970/1261) -30.5 E19-E02 N09-N20
      6 1239 1147 (N11W14 L= 30)   232     72        1 (1208/1119-1263/1168) -45.0 W12-W15 N09-N12
      7 1234 1082 (N07W13 L= 29)   385     68        3 (1204/1055-1271/1107) -14.0 W11-W15 N05-N08
      8 1585  817 (S12W36 L= 52)   447    101        1 (1546/ 791-1633/ 849)  23.2 W33-W40 S13-S10
      9  328  508 (S27E63 L=313)  1056    158        0 ( 286/ 438- 362/ 567)  84.3 E66-E60 S32-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