Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1037 1875 (N58W47 L=328)  1899    248        3 ( 966/1792-1119/1960)   4.3 W31-W85 N53-N60
      2  640 1645 (N55E09 L=272)   569     94        2 ( 603/1617- 687/1677)  11.3 E13-E03 N53-N57
      3 1620 1517 (N15W55 L=335)   372     47        0 (1597/1499-1653/1537) -26.6 W52-W57 N13-N16
      4  604 1426 (N40E17 L=264)   282     52        0 ( 572/1410- 636/1441)  -9.5 E19-E13 N39-N41
      5 1798  973 (S20W49 L=330)   546    173        0 (1766/ 893-1834/1041)  31.0 W44-W54 S23-S16
      6 1577  750 (S26W27 L=307)   265     58        1 (1551/ 728-1605/ 774) -36.9 W24-W28 S27-S23
      7 1396  667 (S25W12 L=293)   614    285        0 (1310/ 565-1504/ 729)  14.0 W04-W20 S28-S22
      8 1264  478 (S33E01 L=279)   574    152        1 (1189/ 462-1327/ 507) -20.6 E05-W03 S35-S29
      9 1612  438 (S50W28 L=308)   733     72        2 (1582/ 400-1635/ 470) -56.3 W24-W31 S55-S46
     10  904  413 (S27E26 L=255)   564    131        0 ( 879/ 361- 940/ 464)  32.0 E29-E21 S29-S23
  

    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