Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  617 1805 (N30E55 L= 31)  1749    277        3 ( 588/1739- 650/1871) -50.2 E60-E50 N25-N37
      2 1605 1441 (N35W26 L=113)  2824    508        4 (1452/1229-1715/1597) -32.6 W11-W37 N24-N42
      3  367 1326 (S05E50 L= 37)  4842    795        6 ( 153/1027- 510/1568)  72.6 E62-E41 S25-N10
      4  280 1525 (N05E70 L= 17)  1139    265        2 ( 246/1500- 320/1559)  -7.4 E78-E63 N03-N07
      5  900 1436 (N13E18 L= 69)   485    136        1 ( 833/1410- 982/1470)  32.0 E22-E14 N10-N16
      6 1289 1258 (N13W08 L= 95)  1777    413        3 (1135/1062-1390/1401) -28.6 E03-W19 N04-N17
      7  540 1089 (S15E30 L= 57)   379    145        1 ( 495/1029- 608/1150)  66.0 E32-E27 S19-S10
      8 1557  662 (S11W43 L=130)   274     85        0 (1544/ 617-1570/ 704) -51.3 W41-W45 S13-S08
      9 1553  450 (S22W54 L=140)   745    148        0 (1496/ 435-1613/ 467)  21.8 W48-W58 S24-S19
  

    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