Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  766 1760 (N31E42 L= 60)   336     48        1 ( 742/1740- 793/1781) -15.9 E45-E38 N29-N31
      2  490 1688 (N19E59 L= 44)  1008    124        1 ( 443/1631- 527/1725) -49.4 E65-E52 N15-N22
      3 1116 1377 (N15E04 L= 98)   515    205        0 (1074/1290-1155/1467) -45.0 E09-W00 N10-N19
      4 1240 1344 (N16W04 L=106)   432    119        0 (1192/1323-1302/1371)  15.9 W00-W07 N14-N18
      5 1342 1254 (N14W12 L=114)   252    102        0 (1325/1204-1365/1304) -45.0 W09-W14 N11-N16
      6 1750  766 (N00W54 L=156)   720    102        1 (1727/ 713-1769/ 818) -23.2 W50-W57 S01-N02
      7 1046  518 (S35W18 L=120)  1810    447        2 ( 896/ 471-1250/ 598)  45.0 W08-W26 S42-S24
      8 1462  499 (S23W45 L=148)   457     59        1 (1430/ 472-1492/ 524)   9.5 W42-W48 S24-S21
      9 1544  341 (S28W65 L=168)   511     88        3 (1525/ 301-1568/ 370) -45.0 W61-W70 S30-S25
     10 1122  296 (S46W38 L=140)   411    109        0 (1066/ 285-1177/ 307)  35.5 W33-W41 S48-S43
  

    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