Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  513 1494 (N31E41 L= 37)   707    105        2 ( 456/1469- 571/1511)  -4.8 E47-E35 N29-N32
      2  933 1361 (N24E07 L= 71)   424    139        0 ( 868/1350-1011/1374)   0.0 E11-E01 N23-N24
      3  701 1261 (N17E22 L= 56)   231     72        0 ( 669/1245- 739/1273)  14.0 E24-E19 N15-N17
      4 1500 1200 (N15W31 L=109)   267    108        1 (1448/1181-1539/1216)   0.0 W27-W34 N13-N15
      5  154  903 (S09E70 L=  8)  1430    108        4 ( 122/ 857- 179/ 956) -45.0 E76-E66 S11-S05
      6  193  821 (S14E65 L= 13)   889    171        4 ( 163/ 765- 230/ 878)  33.7 E71-E59 S17-S09
      7 1424  841 (S08W26 L=104)   290     47        1 (1398/ 812-1450/ 864)  45.0 W24-W27 S09-S06
      8 1662  732 (S14W46 L=124)  1158    201        2 (1601/ 672-1731/ 819)  45.0 W41-W51 S18-S08
      9  384  695 (S20E46 L= 33)  1020    185        1 ( 350/ 613- 436/ 777) -79.7 E48-E42 S25-S14
     10 1077  494 (S31W05 L= 84)   547    142        1 (1056/ 425-1102/ 555) -78.7 W03-W07 S36-S26
     11 1125  374 (S41W10 L= 89)   369     67        1 (1098/ 347-1151/ 406)  53.1 W08-W12 S43-S38
  

    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