Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1697 1426 (N09W55 L=314)   781    174        3 (1669/1366-1725/1492)  54.5 W51-W60 N04-N11
      2 1641 1400 (N09W49 L=308)   287     59        2 (1623/1376-1666/1432)  56.3 W47-W50 N07-N10
      3 1605 1257 (N02W41 L=300)   769    191        3 (1546/1171-1678/1339) -65.6 W39-W44 S03-N08
      4  522  939 (N10E32 L=228)   308     98        0 ( 484/ 914- 557/ 959) -21.8 E34-E29 N07-N11
      5  680  827 (N00E24 L=235)   635    126        1 ( 605/ 801- 744/ 861) -26.6 E28-E20 S01-N03
      6  301  770 (N05E53 L=206)  1192    144        3 ( 254/ 722- 347/ 816) -35.5 E58-E48 N01-N07
      7  451  472 (S15E53 L=207)  1404    331        4 ( 358/ 412- 556/ 523) -39.8 E60-E46 S20-S10
      8  430  409 (S19E58 L=201)   994    141        4 ( 379/ 380- 476/ 440) -39.8 E63-E54 S21-S16
      9 1115  385 (S37E12 L=247)   272     29        0 (1094/ 361-1133/ 404) -90.0 E13-E10 S39-S35
  

    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