Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  722 1837 (N58E52 L=309)  1748    134        2 ( 681/1772- 773/1871) -47.7 E61-E44 N50-N61
      2  857 1608 (N39E19 L=342)  1011    192        1 ( 798/1570- 948/1648) -15.9 E24-E10 N36-N41
      3 1314 1452 (N31W18 L= 18)   288    130        0 (1295/1404-1346/1495) -71.6 W16-W19 N28-N34
      4 1095 1450 (N30W02 L=  2)   290     74        1 (1075/1431-1126/1481) -45.0 W00-W03 N28-N31
      5  233 1418 (N20E72 L=289)  1292     94        4 ( 209/1374- 261/1461) -14.0 E83-E64 N17-N22
      6 1484 1401 (N29W31 L= 31)   297     72        0 (1454/1372-1512/1431) -45.0 W28-W32 N27-N30
      7  613 1344 (N19E30 L=331)   437    161        0 ( 549/1284- 657/1389) -26.6 E35-E25 N15-N21
      8 1706 1180 (N16W47 L= 48)   379     46        1 (1690/1145-1728/1207)  76.0 W46-W49 N13-N17
      9 1139 1086 (N07W07 L=  7)   505    119        1 (1098/1042-1180/1124)  50.2 W04-W09 N03-N09
     10  930  826 (S11E04 L=356)   360    125        0 ( 878/ 789- 986/ 860) -23.2 E07-W00 S12-S08
     11 1702  650 (S17W54 L= 54)   501     44        2 (1678/ 627-1730/ 672)   0.0 W50-W57 S18-S15
     12 1458  559 (S25W35 L= 36)   462    125        0 (1403/ 542-1513/ 590) -11.3 W29-W39 S25-S22
     13 1370  391 (S38W34 L= 35)   306     30        0 (1343/ 380-1397/ 405)   0.0 W31-W36 S38-S36
  

    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