Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  534 1705 (N22E56 L=267)  1506    172        4 ( 484/1643- 575/1784) -70.0 E64-E51 N17-N28
      2  769 1688 (N26E36 L=287)   653     86        2 ( 740/1650- 809/1733) -45.0 E40-E32 N23-N29
      3 1388 1691 (N45W05 L=328)   406     42        0 (1362/1673-1420/1706)  45.0 W02-W07 N43-N47
      4 1184 1392 (N18E00 L=322)   253     51        0 (1161/1364-1206/1417) -33.7 E02-W01 N16-N18
      5 1807 1176 (N24W49 L= 11)   550     76        0 (1778/1130-1827/1209) -31.0 W45-W51 N21-N26
      6  859  653 (S33W00 L=323)   217     74        0 ( 827/ 626- 892/ 674)  45.0 E01-W02 S34-S30
      7 1087  557 (S32W18 L=341)   471     88        1 (1054/ 526-1127/ 599)  -8.1 W14-W21 S33-S29
      8  574  565 (S46E20 L=303)   260     42        0 ( 546/ 548- 598/ 582)   0.0 E22-E17 S47-S45
      9  617  363 (S61E12 L=310)   528     78        1 ( 582/ 339- 651/ 392)   0.0 E17-E07 S63-S57
     10 1518  374 (S28W59 L= 22)   407     68        0 (1483/ 358-1555/ 391)  14.0 W54-W62 S29-S26
     11  786  289 (S61W11 L=333)   588     86        1 ( 750/ 258- 825/ 330)  -5.2 W06-W17 S62-S57
  

    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