Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1116 1888 (N54E34 L=293)   879     81        2 (1084/1863-1156/1915)  71.6 E40-E28 N50-N57
      2  709 1622 (N19E38 L=289)  2591    533        4 ( 663/1382- 758/1858) -49.7 E58-E29 N04-N37
      3 1211 1541 (N29E05 L=322)   405    112        0 (1168/1497-1239/1588) -20.6 E09-E01 N26-N30
      4  930  729 (S25W03 L=330)   619    266        0 ( 881/ 623- 993/ 831) -81.9 W02-W06 S31-S17
      5  581  724 (S36E20 L=307)  1799    332        4 ( 418/ 677- 731/ 776)  30.3 E33-E09 S43-S29
      6 1260  699 (S17W23 L=350)   450    114        0 (1203/ 673-1309/ 722)  31.0 W20-W26 S19-S14
      7 1435  573 (S18W39 L=  6)  1245    280        4 (1368/ 508-1551/ 679)  75.1 W35-W42 S23-S08
      8 1084  607 (S27W17 L=344)   932    190        2 ( 982/ 583-1165/ 640)  28.6 W10-W21 S30-S24
      9  589  350 (S64E11 L=316)  2016    189        9 ( 519/ 314- 667/ 406)  20.9 E23-E02 S68-S58
     10 1486  375 (S27W57 L= 24)   644    143        2 (1418/ 360-1552/ 399)  41.2 W52-W61 S29-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