Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1139 1635 (N41W26 L=338)   269     54        1 (1122/1608-1154/1658)  36.9 W24-W28 N38-N42
      2 1731 1467 (N14W62 L= 14)   930    145        3 (1705/1426-1769/1532)  71.6 W59-W65 N11-N17
      3 1193 1491 (N31W23 L=335)   287     60        1 (1175/1461-1225/1516)   0.0 W20-W25 N28-N32
      4 1783 1408 (N09W65 L= 17)   605    102        1 (1758/1379-1809/1432)   7.1 W61-W69 N07-N10
      5 1503 1213 (N07W33 L=345)   760    214        2 (1482/1096-1522/1312)  67.4 W30-W36 N00-N12
      6  634 1142 (N22E22 L=290)   226     88        0 ( 595/1120- 677/1163)   0.0 E24-E18 N20-N22
      7  140  987 (N19E68 L=244)  1346    289        1 (  94/ 935- 187/1042) -31.0 E72-E62 N15-N21
      8 1607  787 (S21W32 L=344)   717    166        1 (1546/ 705-1653/ 855)  37.9 W26-W36 S24-S17
      9  313  717 (N01E55 L=257)   656     99        1 ( 282/ 675- 348/ 761) -59.0 E57-E52 S01-N04
     10 1371  449 (S37W10 L=322)   309     52        1 (1347/ 427-1395/ 468) -45.0 W08-W11 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