Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  796 1873 (N55E54 L=212)   712     80        1 ( 764/1843- 835/1899)  -9.5 E68-E44 N52-N56
      2 1052 1757 (N50E16 L=250)   457     60        0 (1022/1733-1084/1776)  45.0 E17-E13 N47-N52
      3 1197 1658 (N44W00 L=266)   557    104        0 (1150/1636-1242/1682)   0.0 E03-W04 N42-N46
      4 1677 1377 (N33W44 L=310)   316     15        0 (1655/1363-1699/1392)  21.8 W41-W46 N31-N34
      5 1404 1310 (N23W20 L=286)   242     64        0 (1374/1292-1433/1331) -21.8 W17-W22 N22-N24
      6  180  970 (S17E64 L=202)   630     61        1 ( 158/ 940- 214/ 994)  23.2 E67-E59 S18-S15
      7 1703  677 (S10W54 L=320)   372     85        1 (1675/ 637-1725/ 717)  90.0 W52-W55 S12-S07
      8  903  631 (S26E01 L=265)   300    148        0 ( 836/ 617- 975/ 646)  26.6 E05-W03 S28-S24
      9 1300  475 (S30W30 L=296)   572    124        0 (1237/ 457-1365/ 495)  26.6 W25-W34 S31-S27
  

    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