Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1219 1762 (N55W04 L=167)   504     60        1 (1190/1738-1247/1787)  50.2 W01-W07 N52-N58
      2 1116 1713 (N49E05 L=157)   413    121        0 (1100/1654-1136/1759) -78.7 E06-E01 N43-N53
      3 1502 1712 (N56W40 L=202)   482     39        0 (1483/1685-1527/1731)  26.6 W34-W45 N53-N58
      4  823 1530 (N30E22 L=140)   471    121        3 ( 799/1473- 843/1584) -49.4 E24-E18 N26-N33
      5  372 1466 (N19E56 L=106)  1586    235        6 ( 329/1378- 410/1545) -61.4 E62-E51 N13-N24
      6  294 1440 (N16E64 L= 98)  1213    176        3 ( 259/1374- 334/1517)  90.0 E67-E59 N11-N21
      7 1532 1407 (N31W32 L=194)  1165    300        1 (1418/1325-1647/1475) -20.2 W22-W41 N27-N34
      8  840  939 (S07E10 L=152)   374     95        0 ( 808/ 909- 873/ 965) -26.6 E12-E07 S08-S05
      9  170  903 (S18E67 L= 95)   675     75        0 ( 147/ 869- 192/ 934)  35.5 E71-E62 S20-S15
     10 1240  548 (S26W22 L=184)   328    141        0 (1182/ 517-1296/ 593)  49.4 W18-W24 S29-S22
     11 1700  493 (S24W65 L=227)   609     53        2 (1677/ 467-1719/ 523) -63.4 W61-W68 S25-S21
     12  403  452 (S47E51 L=111)  1176    138        3 ( 365/ 425- 455/ 498) -11.3 E58-E42 S49-S43
  

    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