Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  993 1756 (N48W12 L=164)   338     31        0 ( 974/1729-1016/1776)  36.9 W10-W14 N46-N50
      2  766 1550 (N36E10 L=142)   323     89        0 ( 733/1510- 797/1585)  38.7 E12-E07 N34-N38
      3 1475 1266 (N08W31 L=183)   245     67        0 (1440/1235-1502/1286) -53.1 W29-W32 N05-N09
      4 1690  944 (S14W42 L=194)   275     60        1 (1675/ 911-1702/ 981)  63.4 W41-W43 S15-S11
      5 1394  888 (S13W21 L=173)   230     50        0 (1376/ 855-1413/ 915)  56.3 W19-W22 S14-S11
      6 1322  733 (S21W14 L=166)  1147    382        0 (1155/ 670-1459/ 824)   5.4 W03-W24 S24-S17
      7  426  649 (S14E46 L=106)   702    107        1 ( 387/ 600- 454/ 699)  50.2 E49-E42 S16-S10
      8  878  650 (S20E15 L=137)   287     76        0 ( 856/ 614- 906/ 680)  45.0 E16-E12 S22-S18
      9  225  554 (S17E75 L= 77)  1574    101        3 ( 198/ 496- 245/ 602)  45.0 E80-E70 S20-S13
     10  542  398 (S31E48 L=104)  1917    303        0 ( 453/ 340- 637/ 498) -67.4 E54-E38 S35-S23
  

    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