Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  429 1489 (N31E50 L=179)   345     67        0 ( 407/1467- 457/1522) -26.6 E53-E46 N29-N32
      2  535 1189 (N12E33 L=196)   225     86        0 ( 518/1151- 548/1233) -59.0 E34-E31 N09-N14
      3 1768  952 (S01W53 L=282)   410     68        1 (1744/ 928-1796/ 983) -18.4 W50-W56 S02-N00
      4  428  706 (S19E42 L=187)   314    104        0 ( 381/ 692- 475/ 721)   7.1 E45-E37 S20-S17
      5 1703  653 (S20W52 L=281)   366     49        0 (1684/ 629-1729/ 675)  36.9 W50-W54 S21-S18
      6  430  613 (S25E44 L=185)   534    135        0 ( 413/ 552- 465/ 665)  53.1 E46-E39 S29-S21
      7 1387  583 (S25W26 L=255)   262    113        0 (1352/ 540-1417/ 627)  63.4 W24-W27 S27-S21
      8  804  515 (S31E15 L=214)   275     51        1 ( 771/ 500- 835/ 536)  11.3 E17-E12 S31-S28
      9 1140  368 (S41W11 L=240)   520     43        0 (1111/ 342-1174/ 391)  36.9 W09-W14 S43-S39
     10  887  328 (S46E10 L=219)  1017    193        0 ( 796/ 284- 948/ 364) -24.8 E17-E04 S49-S42
  

    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