Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1405 1392 (N23W31 L=256)   242     78        0 (1364/1359-1430/1426) -68.2 W29-W32 N20-N25
      2  929 1220 (N19E04 L=222)  1812    468        8 ( 708/1143-1116/1309) -32.0 E17-W07 N12-N27
      3 1539 1266 (N13W37 L=262)   360    106        1 (1498/1230-1578/1307) -56.3 W34-W39 N10-N16
      4 1605 1245 (N11W42 L=267)   266     88        1 (1580/1217-1624/1281)  36.9 W39-W43 N09-N12
      5 1390  860 (S09W21 L=246)   471    119        3 (1344/ 836-1450/ 893)   0.0 W17-W24 S10-S06
      6 1417  787 (S14W22 L=247)   945    278        4 (1279/ 679-1516/ 892)  30.5 W11-W29 S18-S08
      7 1864  786 (S23W61 L=286)  1620    270        3 (1807/ 704-1923/ 892)  37.6 W54-W69 S27-S17
      8 1250  786 (S11W11 L=236)   233    106        0 (1215/ 757-1271/ 841)  45.0 W08-W12 S13-S08
      9 1092  775 (S10W01 L=226)   228    110        3 (1067/ 737-1125/ 817) -90.0 W00-W02 S12-S07
     10  296  752 (S04E56 L=170)   867    168        4 ( 246/ 720- 349/ 792) -32.0 E60-E50 S06-S01
     11  395  606 (S14E50 L=175)   797    206        3 ( 375/ 512- 413/ 678)  61.4 E53-E46 S20-S09
  

    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