Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1461 1674 (N47W11 L= 77)   438     47        0 (1430/1655-1494/1694)   8.1 W07-W14 N45-N49
      2 1532 1559 (N40W19 L= 85)   883    109        3 (1486/1512-1569/1607)  -7.1 W14-W22 N36-N44
      3  234 1462 (N01E71 L=355)   692     43        1 ( 210/1438- 264/1479)  -5.2 E77-E66 N00-N02
      4 1009 1425 (N15E11 L= 55)   274    117        0 ( 982/1366-1026/1469) -45.0 E14-E08 N11-N17
      5 1630 1218 (N21W32 L= 98)   271     26        1 (1605/1197-1650/1237)   0.0 W29-W33 N19-N22
      6 1896 1140 (N27W61 L=127)  1483    220        0 (1866/1028-1916/1236) -76.0 W55-W64 N19-N31
      7  217 1136 (S18E60 L=  7)  1107    178        1 ( 155/1096- 293/1170)  14.0 E68-E51 S20-S15
      8 1834 1113 (N22W52 L=119)   500     79        0 (1820/1068-1851/1161) -59.0 W50-W54 N19-N24
      9 1261 1121 (N04W11 L= 77)   233     81        0 (1235/1099-1295/1157) -14.0 W08-W12 N03-N05
     10 1323 1049 (N02W16 L= 82)   244     74        1 (1284/1030-1358/1069)  26.6 W13-W17 N00-N04
     11  179  877 (S34E64 L=  3)   892    104        1 ( 127/ 849- 228/ 913)   0.0 E73-E56 S35-S32
  

    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