Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1022 1405 (N15E09 L=256)   338    112        3 (1005/1353-1048/1446) -80.5 E10-E06 N11-N17
      2 1374 1389 (N22W13 L=277)   455     76        1 (1338/1350-1407/1421) -21.8 W10-W15 N20-N23
      3 1414 1299 (N17W17 L=281)   613    102        2 (1362/1266-1450/1343)  71.6 W14-W19 N14-N20
      4  281 1140 (S14E52 L=212)   508     88        0 ( 246/1111- 312/1173)  38.7 E54-E49 S15-S11
      5 1427 1082 (N05W22 L=286)   224     82        0 (1394/1047-1462/1109) -11.3 W19-W24 N04-N06
      6 1857  902 (N07W59 L=324)   957    186        1 (1834/ 819-1897/ 955) -60.3 W56-W64 N02-N10
      7 1282  739 (S17W22 L=286)   322     73        2 (1239/ 720-1323/ 760)  11.3 W19-W24 S18-S15
      8 1581  413 (S27W56 L=320)  1314    232        3 (1479/ 382-1674/ 450)  34.7 W50-W63 S31-S21
      9  864  206 (S66W21 L=285)   518     86        3 ( 840/ 177- 887/ 230) -29.1 W16-W25 S69-S63
  

    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