Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1723 1541 (N13W67 L=131)  1161     75        4 (1693/1500-1755/1581)  71.6 W62-W71 N09-N15
      2 1543 1535 (N18W49 L=114)   324     70        3 (1510/1512-1571/1555) -63.4 W47-W51 N16-N20
      3 1179 1462 (N25W22 L= 86)  1796    640        5 ( 928/1412-1487/1513) -33.3 W08-W40 N13-N35
      4  289 1475 (N52E56 L=  8)  1739    176        3 ( 204/1440- 361/1507)   0.0 E77-E43 N48-N53
      5  863 1392 (N30W01 L= 65)   238     57        0 ( 838/1371- 895/1413) -71.6 E01-W02 N28-N32
      6  411  999 (N18E39 L= 25)  1334    385        2 ( 378/ 855- 473/1181)  81.5 E42-E34 N09-N29
      7 1832  944 (S23W53 L=117)   432     60        0 (1806/ 913-1853/ 978)  21.8 W49-W55 S24-S21
      8 1260  671 (S21W04 L= 68)   476    242        0 (1140/ 657-1369/ 701) -21.0 E02-W11 S23-S17
      9  446  490 (S12E54 L= 10)   462     74        0 ( 419/ 459- 483/ 517)  14.0 E57-E49 S13-S10
     10  654  353 (S25E46 L= 18)   609    139        0 ( 605/ 310- 681/ 411)  21.0 E53-E40 S26-S21
     11  925  351 (S32E28 L= 36)   502    127        1 ( 901/ 309- 954/ 383)  29.7 E31-E24 S35-S29
  

    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