Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1103 1878 (N58E32 L=219)  1045     83        1 (1076/1835-1127/1920) -40.6 E41-E27 N52-N64
      2  744 1739 (N33E43 L=208)   422     77        0 ( 727/1708- 766/1780) -32.0 E47-E39 N31-N36
      3 1262 1669 (N44E03 L=248)   605     70        1 (1237/1631-1290/1715) -39.8 E06-W00 N40-N47
      4 1577 1589 (N49W28 L=279)   777     88        2 (1533/1562-1619/1612)   0.0 W22-W33 N46-N51
      5 1700 1473 (N44W41 L=292)  1138    212        2 (1646/1403-1762/1549) -41.6 W35-W49 N40-N48
      6 1613 1262 (N26W32 L=283)  1832    531        2 (1470/1046-1756/1423) -28.9 W18-W47 N16-N32
      7  173 1006 (S24E65 L=186)  1739    201        3 ( 138/ 905- 213/1082)  67.4 E71-E58 S30-S18
      8 1085  619 (S25W15 L=266)  1354    287        2 ( 963/ 570-1225/ 661)  40.6 W08-W22 S31-S19
      9 1306  533 (S24W32 L=283)   318    119        0 (1260/ 491-1342/ 575)  74.1 W30-W32 S27-S20
     10  801  526 (S39E01 L=250)   250     72        1 ( 776/ 509- 838/ 553) -14.0 E02-W02 S39-S37
     11 1450  350 (S31W55 L=306)  1896    219        3 (1338/ 326-1539/ 384)  14.0 W44-W64 S33-S27
     12  972  159 (S63W44 L=295)  2053    175        2 ( 937/ 100-1029/ 228)  -9.9 W29-W75 S67-S56
  

    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